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A systematic method of summing the corrections to the renormalon residue arising from higher 
order renormalons is discussed. 



1 Introduction 



It is well known that the weak coupling expansion in field theory is not convergent but 
asymptotic. Its coefficients diverge factorially in the order of perturbation. One source 
of the divergence is renormalon. For example, when a Green's function such as the Adler 
function D(a) is expanded as 

D(a) =Y j a n a n+ \ (1) 
the infrared renormalon causes the coefficients to grow as 



Kn\n- 2/3l/ ffi (-2//? )" n , for n -► oo, (2) 



where /3q,P\ are the first two coefficients of the (3 function. 

The constant K is an all order quantity: it gets contribution not only from the well- 
known renormalon diagram, a chain of one loop bubbles, but also from infinitely many higher 
order renormalon diagrams. The purpose of this talk is to show summing these corrections 
is manageable, and that K can be expressed in a calculable, convergent sequence. 

2 Renormalons in QED 

For our purpose it is convenient to consider the Adler function defined by 



oo 

= const. + a n (fj, 2 /q 2 )a(f,) n+1 (3) 



n=0 



where II is the vacuum polarization function of the electromagnetic current and \i is the 
renormalization scale. 

The large order behavior in Eq. |2| arises from an exchange of the all-order Gell-Mann- 
Low (GL) effective charge S defined byi 

a « a '"»' t » = 1 + ntw, t ) (4) 



"Presented at Vllth Blois Workshop, Seoul, Korea, June 10-14, 1997. 



1 



where t = k 2 /fi 2 . Therefore, as far as we are concerned with the infrared renormalon 
associated with the asymptotic behavior, we can write the Adler function as 



D(a( t i),fi 2 /q z )= f(k') a (k 2 )dk z (5) 
Jo 

where a.(k 2 ) denotes the GL effective charge and 

f{k2) = fOTk2 ~* (6) 

with ef denoting the electric charge. 

The divergence in Eq. ^ causes a singularity in Borel plane. The Borel transform of 
D(a(fi), /j? /q 2 ) is defined by 

D(a(p) } fJ?/q 2 ) = J™ exp (~^y) 5(6) d 6. (7) 
With the perturbative series in Eq. [l] the Borel transform becomes 

oo 

D(b) = £ ^6" (8) 



n=0 



Substitution of the asymptotic form in Eq. |2| into Eq. |s| gives 

bib) - (9) 

in the neighborhood of the singularity at b = —2//3q. Thus, to calculate K we need to 
evaluate the residue of the singularity. 

Let us now consider the renormalization group equation for the GL effective charge: 



/^ttt + P(a)— a(a(/x), t) = 0. (10) 



Solving the equation we may write the GL effective charge as 

where C(a) is a scheme independent function. The effective coupling A is defined by 

A(a(u),t) = 7 -\ r , (12) 

where p% is the constant term in the one loop vacuum polarization function. The function 
C(a(/c 2 )) may be expanded in an asymptotic series as 

oo 

C( a (k 2 )) = ci ln(a(A; 2 )) + £ a K*; 2 )]*- 1 . (13) 

2 
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The Ci can be easily determined in terms of the coefficients of the (3 function and the vacuum 
polarization function. 

Now define an N th order GL effective charge: 



aW(fc 2 ) 



1 



and introduce a modified Borel transform 

b 







cxp 



A(t) 



iW(b)db. 



(14) 



(15) 



Substituting Eq. 15 into Eq. ||, we can write D^ N \a(fi), fi 2 /q 2 ), which is defined by replacing 
a(k 2 ) in Eq. | with a^/c 2 ), as 



D^ia^l^/q 2 ) 



exp 



bp 



a(p) da 



/3(a) 



{e- bpi /(&)a (7V) (6)}d6, 



where 



f(b) = / f(t) exp (6ft In t) dt. 



(16) 



(17) 



The first IR renormalon singularity arises from the IR divergence in the integral in Eq. [I?]. 
Substituting Eq. ||| into Eq. 17 



f(b) 



A I 



f{t)e bfi ° hlt dt 

<-M 

te bf3 ° lnt dt 

1 



e 2 ft JL A rM 



8ttV Jo 



8vrV 2 + bp 



l + (2 + 6/3 )lnM + 



(18) 



where M is an arbitrary UV cutoff. Notice that the leading renormalon singularity is cutoff 
independent. 

In the renormalization scheme in which the (3 function is given by 



/3(a) 



A)« 2 

1 — Aa 



(19) 



with A = /3i//3o h the Eq. 16 becomes the modified Borel transform introduced by Brown, 
Yaffe, and ZhaiS. Then using the relation between the ordinary Borel transform and the 
modified one, we have the ordinary Borel transform 



=2 ,,4 



-2A//3 )! 



r>( N Mh\ — > -b oPl -Xb Q lnb Q ~(N)( h \ 

{0) Wtt^ 6 [0 °\ 1 + i b p y-2x/fs 



(20) 



where bo = — 2//3q. 



Using Eq. 14 the Borel transform of the effective charge can be obtained without diffi- 



culty. The final expression is given byu: 



ifici r oo / h \k HN-i) , N 



k=0 '"' l=k y i=0 



oo k(JM-l) N , , fe 

EV^ \ " l~ -U n NklCj ,k+l+i 

f-f ^ /t!r(/ + z + l-6ci) 1 j 



fc=0 «=fc i=0 



where 



and 



1 for z = 

ci for i = 1 (22) 
'i — 1) Cj for i > 2, 

yiN-l ni 

h mi = k\J2 ^-1 i+ l ( 23 ) 

with the set {n^} of nonnegative integers satisfying 

AT-i JV-l 

5^ ni • i = Z, 51 ni = /c. (24) 

i=l i=i 

To find the renormaion residue of D^ N \b), we have to evaluate b bci a^ N \b) at the first IR 
renormaion position, bo = —2/(5q. If we directly substitute b in Eq. 21 with bo, the resulting 
large order behavior does not have a finite limit for N — > ocM The reason for this is that 5(6) 
is singular at the UV and IR renormaion positions, and its radius of convergence when it is 



expanded as in Eq. 21 is given by the position at b = 1/(3q of the first UV renormaion, which 



is the closest renormaion to the origin in the Borel plane. Therefore we cannot substitute 



b with bo in Eq. 21 to correctly evaluate the Borel transform at the first IR renormaion. 

This problem can be avoided by introducing an analytic transform of the Borel plane 
so thaLthe closest renormaion to the origin in the new complex plane is the first IR renor- 
maion H. Because the singularity of a(6) at the IR renormaion is such that it is finite but 
has a divergent derivative 13, we can then express the residue as a convergent series. 

For this purpose, we can take any analytic transform that puts the IR renormaion as 
the closest singularity to the origin. Here we take: 

l - /?oo 

In the z — plane, the closest singularity to the origin is the first IR renormaion at 

zo = -~, (26) 
and all the UV renormalons are pushed beyond z = — 1 on the real axis. 
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Table 1: The first four elements of the sequence for the first IR renormalon residue in QED 





Nf = l 


N f = 2 


JV> = 3 


N f = A 


JV) = 5 


Nf = 100 


Kq 


1 


1 


1 


1 


1 


1 


«1 


1.63 


1.32 


1.21 


1.16 


1.13 


1.00 


K2 


0.71 


1.31 


1.31 


1.27 


1.24 


1.02 


K3 


-1.53 


1.25 


1.41 


1.39 


1.34 


1.02 



To find b bci &( N \b) at the first IR renormalon, we have to substitute b in Eq. 21 with 
the inverse of Eq. ^5] and expand it in Taylor series at z = to order N, and evaluate it at 
z = zq . Then the Borel transform of GL effective charge at the first IR renormalon is given 
by 

N 

EM 



K N 



b bC ^ N \b) 



b=b 



(27) 



Af=0 



where it is straightforward to express qM in terms of Cj. 
Finally, the coefficient K in Eq. |2| is given by 



where 



K N 



K = lim Kn 



e 2 /i 4 

/P T e- b0Pl K N . 



(28) 



(29) 



16vrY 

The first four elements of the sequence is given in Table 1 for different number of fermion 
flavors. 



3 Renormalons in QCD 

In QCD, there is unfortunately no satisfactory definition of renormalization scheme and 
scale invariant effective charge that may be used in the diagrammatic study of renormalon. 
However, if we are only interested in the calculation of the residue, the definition of the 
effective charge is not required. Indeed the calculation is cunningly simple; it only requires 
the strength of the renormalon singularity and the perturbative calculation of D(a). For this 
reason, the method described in the following can be equally well applied to the calculation 
of the UV renormalon residue. 

Consider the Borel transform of the current correlation function in QCD. The renor- 
malon singularity of D(b) in QCD 

6 < 6 > " (!-»/£)■+*■ (30) 

gives the large order behavior 



To calculate the residue D, consider a function 

R(b) = D(b)(l-b/b ) 1+Xb °. 



(32) 



Then because of Eq. 30, we have 



D = R(b ). (33) 

To avoid the first UV renormalon, we introduce a new variable z, as we did in QED, which 
is defined by 



(3 b 



l-Pob 



with its inverse 



b=± 



In the z— plane, the IR renormalon at 



(34) 
(35) 

(36) 



is the closest singularity to the origin, and so the radius of convergence of the Taylor series 
of D(b(z)) at z = is given by the first IR renormalon. 
Now D can be expressed in a convergent series form 



D 



D(b)(l-b/b ) l+Xb ° 



b=b 



£^ [&(*)]" (l-6(.)/M 1+A6 ° 

\n=0 H - J 
oo 

n=0 



z=z 



(37) 



where it is straightforward to find r n in terms of the perturbative coefficients a n . Note that 
the series is convergent even if R(b(z)) is not analytic at z = Zq, because then the radius of 
convergence of the series is given by z = Zq, and R(b(zo)) is finite. 

Using the perturbative calculation of the current correlation function, and D(a), to 
three loop!'!, we have 



R(b(z)) 



16 7T 3 



1.333 - 0.748 z - 0.311 z 2 + 0(z 3 



(38) 



for Nf = 3. This is in the renormalization scheme in which the one-loop renormalization 
point is same as that of MS scheme, and the (3 function is given in the form in Eq. [l^. 
Evaluating this series at the renormalon position at z = zq we have 

1.333 



K 2 
K 3 



0.946 



(A&o)! 

(1.323 -0 748,0) =Q592 
(A&o)! 

(1.323 - 0.748 z - 0.311 zg) 
(A6oT! 



0.494 



(39) 
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Table 2: The first three elements of the sequence for the large order behavior in QCD 





N f = l 


N f = 2 


Nf = 3 


Nf = 4 


N f = 5 




.881 


.904 


.946 


1.018 


1.132 


K 2 


.521 


.546 


.592 


.674 


.813 


K 3 


.592 


.549 


.494 


.411 


.307 



For several other flavor numbers we give K n in Table 2. 
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